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Abstract. We describe a new algorithm for calculating the topological 
degree deg (/, B, 0) where B C R n is a product of closed real intervals 
and f : B —> K" is a real- valued continuous function given in the form of 
arithmetical expressions. The algorithm cleanly separates numerical from 
combinatorial computation. Based on this, the numerical part provably 
computes only the information that is strictly necessary for the following 
combinatorial part, and the combinatorial part may optimize its com- 
putation based on the numerical information computed before. We also 
present computational experiments based on an implementation of the 
algorithm. Also, in contrast to previous work, the algorithm does not as- 
sume knowledge of a Lipschitz constant of the function /, and works for 
arbitrary continuous functions for which some notion of interval arith- 
metic can be defined. 



1 Introduction 

The notion of topological degree was introduced by Jan Brouwer [5] and was 
motivated by questions in differential topology |19I12] . The degree of a contin- 
uous function is an integer, describing some topological properties of it. Degree 
theory has many applications, including geometry [26], nonlinear differential 
equations [171617118] , dynamical systems [13] , verification theory [16] , fixed point 
theory [1] and many others. 

The presented algorithm is able to calculate the degree deg (/, B, 0) of any 
real-valued continuous function defined on a box B such that ^ f(dB) and 
/ given in the form of arithmetical expressions containing function symbols for 
which interval enclosures can be computed [20125) . 

Experimental examples show that low-dimensional examples of simple func- 
tions (up to dimension 10) terminate in reasonable time0 In addition to effi- 
ciency, the algorithm has several advantages over previous work that we now 
describe in more detail. 

Erdelsky proposed in 1973 a method for degree computation, assuming that 
the function in question is Lipschitz with a known Lipschitz constant [5] . Another 
method involved the approximation of a multidimensional integral of a function 
derived from / and its partial derivatives. Thomas Neil published an algorithm 



1 The program is accessible on topdeg.sourceforge.net. 



based on integration in 1975, although the error analysis uses only probabilistic 
methods 23 . Other authors constructed algorithms that cover the boundary dB 
with a large set of (n — l)-simplices and use the information about the signs of 
fj on the vertices of these simplices to calculate the degree in a combinatorial 
way. However, the calculated result is proved to be correct only if a parameter is 
chosen to be sufficiently large [14127128] . Boult and Sikorski developed a different 
method for degree calculation in eighties, but their algorithm also requires the 
knowledge of a Lipschitz constant for / [2J. Later, many algorithms arose where 
the degree was calculated recursively from partial information about / on the 
boundary dB. One of such identities is 

deg(/, J B,0)=deg(/. 1 ,C/,0), 

where = (/2, ...,/„), U is a [d — l)-dimensional open neighborhood of 
{x € dB\ f-,i(x) = fi(x) > 0} in dB. Aberth described an algorithm using 
this formula, based on interval arithmetic pQ, although it was not implemented 
and is rather a recipe then a precise algorithm. Later, Murashige published a 
methods for calculating degree that uses concepts from computational homology 
theory [21]. 

Although a broad range of ideas and methods has been implemented, the 
effectivity of these algorithms decreases fast with the dimension of B. For ex- 
ample, in the Murashige homological method, computation of the degree of 
the identity function f(x) = x takes more then 100 second already in dimen- 
sion 5 [21] Figure 3] . Other approaches were developed that calculate degree of 
high-dimensional examples quickly, provided the functions are of some special 
type. For instance, there exist effective degree algorithm for complex functions 
f:C n ^C n [STT5TTB] . 

We formalized, extended and implemented the rough ideas of Oliver Aberth 
PQ. The algorithm presented here works in a very general setting. It assumes 
that the function / is real valued and continuous, defined in terms of inter- 
val arithmetic. We don't require the function to be differentiable and not even 
Lipschitz. 

From the algorithmic point of view, the main step consists of covering the 
boundary of a d-dimensional set SI with (d — l)-dimensional regions D\, . . . , D m 
where a particular component /; of / has constant sign and then recursively 
gathering the information about the signs of the remaining components of / on 
dDj. All the sets are represented as lists of oriented boxes. They do not have 
to represent manifolds and we allow the boundary of these sets to be compli- 
cated (see Def. U]). ^ n this setting, it is computationally nontrivial to identify 
the boundary dDj of a d-dimensional set embedded in K™ and to decompose 
the boundary to a sum of "nice" sets. Instead of doing this, we calculate an 
over-approximation of dDj that is algorithmically simpler and then prove that 
it has no impact on the correctness of the result. This involves some theoreti- 
cal difficulties whose solution necessitates the development of several technical 
results. 

The second section contains the main definitions needed from topological 
degree theory. Theorem [T] is a fundamental ingredient of our algorithm. Section 



3 describes the algorithm itself and its connection to Theorem [TJ In Section 4, 
we present some experimental results. The last Section contains the proof of 
two auxiliary Lemmas that we need throughout the paper. These proofs do not 
involve deep ideas but are quite long and technical, so we moved them to the 
end of the paper in order to keep the main parts of the text fluent and readable. 

2 Mathematical Background 
2.1 Definitions and Notation 

In this section, we first summarize the definition and main characteristics of the 
topological degree on which there exists a wide range of literature |9I24[ e.g.,]. 
Degree theory works with continuous maps between oriented manifolds, and in 
order to represent these topological objects on computers we will then introduce 
Definitions Q] to [6] Finally, original Theorem Q] will be the main ingredient of our 
algorithm for computing the topological degree. 

Let fl C R™ be open and bounded, / : Q —> R" continuous and smooth (i.e., 
infinitely often diffcrcntiable) in !7, p ^ f(dfl). For regular values p £ R™ (i.e., 
values p such that for all y £ det/'(y) ^ 0), the degree deg (/, p) is 

defined to be 

deg(f,Q,p):= signdet/'(y). (i) 

v&f-Hp) 

This definition can be extended for non-regular values p in a unique way, so 
that a(p) :— deg (/, fi,p) is locally constant on the connected components of 
R n \f(d(2) Ql. 

Here we give an alternative, axiomatic definition, that determines the de- 
gree uniquely. In this approach deg (/, fi,p) is the unique integer satisfying the 
following properties [9124110] : 

1. For the identity function /, deg (/, f2,p) = 1 iff p is in the interior of fl. 

2. If deg (/, fl 7 p) 7^ then f(x) = p has a solution in fl. 

3. If there is a continuous function (a "homotopy") h : [0, 1] X Q — > R" such 
that p i /([0, 1] x dfl), then deg (/, Q,p) = deg (g, Q lP ). 

4. If Q 1 r\Q 2 =0andp^ f{Q\ Utf 2 )), then deg (/, Q,p) = deg(/,/2 1 ,p) + 
deg(/, rj 2 ,p). 

5. The function a(p) := deg (/, fl,p) is constant on any connected component 

of R n \f (on). 

This can be extended to the case where M and N are any compact oriented 
manifolds of the same dimension. Formula ([T|) is still meaningful, f'(y) denotes 
the matrix of partial derivatives of some coordinate representation of /. For 
example, if / is a function from an oriented curve c (i.e., an oriented set of 
dimension 1) to R and / / on the endpoints of c, then deg (/, c, 0) is well 
defined. 

If / : M — > N is a function between two oriented manifolds without boundary, 
then the degree deg (/) is defined to be deg (/, M,p) for any p 6 f{M). The 



axiomatic definition of degree is meaningful for all continuous functions between 
oriented manifolds of the same dimension. A simple consequence of these pro- 
perties is that if p £ f(dfi), then deg (/, H,p) — deg (/ — p, J?, 0), so we will be 
only interested in calculating deg (/, fl 1 0). 

We will represent geometric objects like manifolds, orientation, boundaries 
and functions in a combinatorial way, using the following definitions. 

Definition 1. A k- dimensional box (simply k-box) in K" is the product of k 
non- degenerate closed intervals and n — k degenerate intervals (one-point sets). 
A sub-box of a k-box A is any k box B s.t. B C A. 

Definition 2. The orientation of a k-box is a number from the set {1, —1}. An 
oriented box is a pair (B, s) where B is a box and s its orientation. We say that 
B\ is an oriented sub-box of an oriented box B, if B\ C B, the dimensions of 
B and B\ are equal and the orientations are equal. 

Definition 3. Let B = I\ X I2 X . . . X 7„ be an oriented d-box in W 1 with ori- 
entation o. Let, for every i G {1, . . . , n}, [oj, 6j] = Li. Assume that the intervals 
Lj t , . . . Lj d are non-degenerate, j% < j% < • • • < jd, the other intervals are degen- 
erate (one-point) intervals. For i S {1, . . . , d}, the (d — 1) -dimensional boxes 

Ff := {{xi, . . . , x n ) G B I x n — a ji } and F+ := {(xi, . . . , x n ) <E B \ x n — b ji } 

are called faces of B. Any sub-box of a face is called a sub-face of B. If we choose 
the orientation of F^ to be (— l) i+1 o and the orientation of F^ to be {— l) l o, 
then we call F^ 1 oriented faces of B. An oriented sub-box of an oriented face is 
called oriented sub-face. The orientation of the oriented faces and sub-faces is 
called the induced orientation from the orientation of B . 

Definition 4. An oriented cubical set Q is a finite set of oriented boxes B\, . . . , B^ 
of the same dimension d such that the following conditions are satisfied: 

1. For each i ^ j , the dimension of Bi n Bj is at most (d— 1). 

2. Whenever Bi n Bj = Bij is a (d — 1)- dimensional box, then the orientations 
of Bi and Bj are compatible. This means that Bij has an opposite induced 
orientation as a sub- face of Bi as the orientation induced from Bj . 

The dimension of an oriented cubical set is the dimension of any box it contains. 
If fi is an oriented cubical set, we denote by \ fi\ the set it represents (the union 
of all the oriented boxes contained in fl). 

An oriented cubical set is sketched in Figure [T] An immediate consequence 
of the definition is that each sub-face F of a box B in an oriented cubical set 
fl is a boundary sub-face of at most two boxes in £1. Note that an oriented 
cubical set does not have to represent a manifold, because some boxes may have 
lower-dimensional intersection, like B\ and -B4 in Figure [TJ 

Definition 5. An oriented boundary of an oriented d-dimensional cubical set 
Q is any set of (d — I) -dimensional oriented boxes df2, such that 



1. Any two boxes in dfl have intersection of dimension at most d — 2. 

2. For each Fq S dfl, and each [d — 1)- dimensional sub-box F' of Fg, there 
exists exactly one box B £ fl such that F' is an oriented sub- face of B. 

3. dfl is maximal, that is, no further box can be added to dfl such that condi- 
tions]]] and\^ still hold. 

An oriented cubical set and its oriented boundary are sketched in Figure [5J 
Geometrically, this definition describes the topological boundary of an oriented 
cubical set fl and we denote the union of all oriented boxes in dfl by \dfi\. 
Clearly, d\fl\ — \dfl\, the meaning of the left hand side being the topological 
boundary of the set \fl\. Note that if fl is a <i-dimensional oriented cubical set 
and dfl an oriented boundary of fl, then each sub-face x of some box in fl s.t. 
x (~1 \dfl\ is at most (d — 2)-dimensional, is a sub-face of exactly two boxes in fl 
with opposite induced orientation (see B\i in Figure [4]). 

An oriented boundary of an oriented cubical set does not have to form an 
oriented cubical set, because the second condition of Definition @] may be violated 
(for a counter-example, see Fig. [3] where the 1-boxes a and c have O-dimensional 
intersection but not compatible orientations). 

The notion of topological degree can be naturally generalized to oriented 
cubical sets. So, if / is a continuous function from a d-dimensional oriented 
cubical set fl to R d such that ^ f(d\fl\), then deg (/, fl, 0) is well-defined, 
extending the definition of deg (/, fl, 0) for oriented manifolds fl. 

Finally, we will represent functions as algorithms that can calculate an over- 
approximations of f{B) for any given box B. 

Definition 6. Let fl C R™. We call a function f : Q — > M interval-computable 
if there exists an algorithm that, for a given box B C fl with rational endpoints 
computes a closed (possibly degenerate) interval I(f)(B) such that 

- I(f)(B) 2 {f(x) I x E B}, and 

— for every e > there is a 5 > such that for every box B with diam(B) < 5, 
for all y G I(f)(B), there is an x G B such that \f(x) — y\ < e. 

We call a function f — (/1, . . . , /„) : fl — > R™ interval- computable iff each fi is 
interval-computable. In this case, the algorithm 1(f) returns a tuple of intervals, 
one for each fi . 

Usually such functions are written in terms of symbolic expressions contain- 
ing symbols denoting certain basic functions such as rational constants, addition, 
multiplication, exponentiation, and sine. In that case, the first property above 
follows from the fundamental theorem of interval arithmetic, and the second 
property from Lipschitz continuity of interval arithmetic (e.g., Theorem 2.1.1 in 
Ncumaier's book [22 ). We will use interval computable functions and expres- 
sions denoting them interchangeably and assume that for an expression denoting 
a function /, a corresponding algorithm /(/) is given. 

2 For an oriented cubical set fl, one can define an oriented manifold O e := {x £ 
int(\fl\) I dist(x, dfl) > e} for a small enough e and define the degree to be 

de g (/,rr,o) 
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Fig. 2. Oriented cubical set Q — {Bi,B2} with orientation of both boxes 1 and its 
oriented boundary dO — {ai, . . . , an} with orientation indicated by the arrows. For 
example, ai has orientation —1 (the arrow goes in the opposite direction as the vertical 
axis), a,2 has orientation 1 etc. 



2.2 Main Theorem 



Now we define the cubical information we use to compute the degree, and prove 
that it is both necessary and sufficient for determining the degree. 

Definition 7. An d- dimensional sign vector is a vector from { — ,0, +} d . 

Let S be a set of oriented (d— l)-boxes. A sign covering of S is an assignment 
of a d-dimensional sign vector to each a £ S . For a sign covering SV and a G S 
we will denote this sign vector by SV a , and its i-th component by {SV a )i- 

A sign covering is sufficient if each sign vector contains at least one non-zero 
element. 

A sign covering is a sign covering wrt. a function / : (U ae s a) — > M. d 
with components (/l, • • • , /d), if f or every oriented box a G S and for every 
i G {1, . . . , d}, {SV a )i 7^ implies that fi has constant sign (SV a )i on a. 

In the following we will often recursively reduce proofs/algorithms for d- 
dimension oriented cubical sets, to proofs/algorithms on their oriented boundary. 
Since — as we have already seen — an oriented boundary of an oriented cubical 
set does not necessarily have to form an oriented cubical set, we will need the 
following lemma that will allow us to decompose this oriented boundary again 
into oriented cubical sets: 

Lemma 1. Let Q be an d-dimensional oriented cubical set, df2 an oriented 
boundary of fi, SV a sufficient sign-covering of dQ with respect to f : \ fl\ — > M. d 
and assume that for each a G dfl, SV a has exactly one nonzero component. 
Let A v<sl := {a G df2\(SV a )i> = s'} for each I' G {l,...,d} and s' G {+,-}. 
Then there exist oriented cubical sets Di, . . . , D m and corresponding oriented 
boundaries dD\, . . . , dD m s.t. the following conditions are satisfied: 

1. Uj E{1 ,... jm} .Dj = dfl, 

2. A n Dj = for i f j, 

3. For each i, there exists l(i),s(i) such that Di C Aiu\ s u\, 

4- Each b G dDi is a sub-face of some a G Ai^ s i where I' ^ l(i). 

The lemma is illustrated in Figure [3] The proof of this Lemma is technical and 
we postpone it to the Appendix in order to keep the text fluent. 

Theorem 1. Let Q be an oriented d-dimensional cubical set, dft an oriented 
boundary of SI and f : \ f}\ — > M. d a continuous function with components (/i, . . . , fd) 
such that ^ f(\dfl\). Then a sign covering SV of dQ wrt. f determines the 
degree deg (/, fl, 0) uniquely if and only if it is sufficient. 

Proof. We first prove that sufficiency of the sign covering implies a unique degree. 
We proceed by induction on the dimension of Q. If fl is 1-dimensional oriented 
cubical set ao, then deg (/, f?,0) = ^ (sign (/(&)) — sign (/(a))) is determined 
by the sufficient sign covering of df2 wrt. /. Let d > 1. For each box a G dSl, 
choose an index i(a) such that (S'V a ) i ( a ) =: s a ^ 0. For all I' G {1, ...,d} 
and s' G {+,—}, let Ai>, s > := {a G df2\i(a) — V , s a — s'}. It follows from 



Lemma [T] that we may decompose dQ into oriented cubical sets Dj and oriented 
boundaries dDj, j — 1, . . . , m such that Di C Aiu\ s u-\ for unique l(i),s(i) and 
each x £ dDi is a sub-face of some & e ^i', s ' where /' For each 

define f^ v := (f 1} . . . , fv-\,fv+u • - • > U)- Then (/ f^(\dDi\) and the degree 
deg(/-,i(i),-Dj,0) is defined. Let I € and s 6 {+, — } be arbitrary. It 

follows from [TSJ Theorem 2.2] and [371 Section 4.2] that 

deg(/,r2,0 d ) = S (-l) i+1 £ deg(/^,A,O d -i) (2) 

i; and s(i)—s 

where 0^ G K. fc is the fc-dimensional zero. 

For each set Di from the sum on the right hand side, fnj\ has sign s(i) on 
Di. Each x £ dDi is a sub-box of some b G vl{/ )S < where /' so we may 

assign a new sign vector for x by deleting the Z(i)'th component from S'Vj,. In 
this way, we define a sufficient sign covering of dDi wrt. and the degree 

deg (/, B, 0) can be then calculated recursively using 

Now assume that the sign covering of dQ is not sufficient. We will prove that 
in this case, the degree is not uniquely determined. 

Let F € dQ be a (d— l)-dimensional box such that SVf — (0, . . . , 0). Choose 
mEZto be arbitrary. We will construct a function G : — > R d such that the 
sign covering of dQ is a sign covering with respect to G and deg (G, 0) = m. 

Denote the oriented manifold with boundary dQ \ F° by Si . dQ is a union of 
the oriented manifolds Si and F, the boundaries dF and 95*1 are equal with op- 
posite orientations, homeomorphic to the sphere S d ~ 2 . The degree deg (/, Q, 0) = 
deg (/) where / = f/\f\ : dQ ->• S^ 1 C R d is a map to the sphere. Let p e S^ 1 
be such that p f(dSi), let a = deg (/, Si,p) and m' = m — a. We construct a 
map g : F —> S^ 1 such that deg (g, F,p) — m' . The homotopy group TTk(S l ) = 
for k < I, so each map from a (d — 2)-sphere to the (d — l)-sphere is nomotopic 
to a constant map. Let us define g\ = f on dF ~ S d ~ 2 . Then gi : dF — > S"* -1 is 
homotopic to a constant map. There exists a sub-box F' <Z F and a continuous 
extension g 2 : -F \ (-P') — ^ of gi such that 52 = 9% = f on and 32 is 

constant on <9F' ~ S" 1-2 . Using the fact that 7r c ;_ 1 (5 d_1 ) = Z, there exists a 
map h : S 4 ' 1 S^" 1 of degree m'. It follows from the identity S"* -1 ~ F' /dF' 
that we can extend g 2 to a map 53 : F — > S"* -1 such that deg(g 3 ,F,p) = m'. 
Finally, extend 53 to a map g : dQ — > S^ -1 by g = f on Si. Then 

deg ( 5 ) = deg (g, Si,p) + deg ( 3 , F,p) = a + m' = m. 

Let i : R d be the inclusion. Multiplying i o g by some scalar valued 

function, we can obtain a function g' : dQ — > M d such that g' = f on Si. 
Extending g' : dQ — > R d to a continuous G : J? — > R d arbitrarily (this is 
possible due to Tietze's Extension Theorem), we obtain a function G such that 
the original sign covering is a sign covering of dQ wrt. G and deg (G, X?, 0) = m. 
This completes the proof. □ 



3 Algorithm description 



3.1 Informal Description of the Algorithm 

We describe now our algorithm for degree computation of an interval computable 
function. If / : B — > R n is an interval computable function nowhere zero on 
the boundary dB, then the corresponding interval computation algorithm /(/) 
from Definition [5] may be used to construct a sufficient sign covering of dB 
wrt. /. This sign covering will be represented as a list of oriented boxes and 
sign vectors. The main ingredient of the algorithm is Equation ^ from the 
proof of Theorem Q] For some index I and sign s, we select all the boxes a 
with (SV a )i — s. From LemmaQ] we know that these boxes form some oriented 
cubical sets Di, . . . , D m . Then a new list of oriented boxes and sign vectors 
is constructed that cover the boundaries dDj of Dj. The new sign vectors are 
assigned to these boundary boxes in such a way that the new list is sufficient 
and wrt. /-,; = (fx, . . . , fi+i, • • • , f n ). Equation ^ is used for the recursion 
step. 

We work with list of oriented boxes and sign vectors rather then with sets, 
because it will be convenient for our implementation to allow an oriented box 
to be contained in the list multiple times. However, we will usually ignore the 
order of the list elements (i.e., the algorithm actually is based on multi-sets 
which we implement by lists). For two lists Lx and L2, we denote by L\ + L2 
the concatenation of L\ and L 2 . We will also use the corresponding symbol 
for list concatenation. We use the notation a G L if a is contained in L at least 
once. If Lx is a sub-list of L, we denote by L \ Lx the list L with the sub-list Lx 
omitted. 

Now we define a version of the notion of sign covering based on lists: 
Definition 8. A sign list (of dimension d) is a list of pairs consisting of 

— an oriented d-box, and 

— a corresponding [d+ 1)- dimensional sign vector. 

A sign list is sufficient, iff each sign vector contains at least one non-zero ele- 
ment. A sign list of dimension d is a sign list wrt. a function / : [J aeL a — > M. d+1 
iff for each element a 6 L and corresponding sign vector SV a = (si, . . . , Sd+i), 
for all i G {1, . . . , d + 1}, Sj ^ implies that fi has sign Si on a. 

By misuse of notation, we will sometimes refer to the elements of a sign list 
as pairs consisting of an oriented box and a sign vector, and sometimes just as 
an oriented box. 

The basic ingredient of the algorithm is a recursive function Deg with input 
a sufficient sign list and output an integer. This function involves no interval 
arithmetic and is purely combinatorial. If the input is a sufficient sign list L wrt. 
/ such that the boxes in L form an oriented boundary of an oriented cubical 
set fi, then this function returns deg(f, fi, 0). If the Deg function input is a 0- 
dimensional sign list L, then the output J2 p ^l orientat ' l ° n ^ XSp is returned. This 



is compatible with the the formula for degree of a function on an oriented edge, 
deg (/, at, 0) = sign/(fc) " sign/(a) . 

If the input consists of oriented <i-boxes and sign vectors of dimension d + 1 
for d > 0, we choose I G {1, 2, . . . , d + 1} and s 6 {+, — } and compute a list 
of boxes L sel (the selected boxes) having s as the I'th component of the sign 
vector. We split the boundary faces of all selected boxes until each face £ of a 
selected box a is either contained in some non-selected box or the intersection 
of x with each non-selected box is at most (d — 2)-dimensional. For each face x 
of a selected box a that is sub-face of some non-selected box b, we delete the l-th 
entry from the sign vector of b and assign this as a new sign vector to x. The 
list of all such oriented (d— l)-boxes and their sign vectors is denoted by faces. 
This is a sufficient sign list wrt. /-,; and s (— Deg(faces) is returned. 

The choice of / and s has no impact on the correctness of the algorithm but 
can optimize its speed. We choose I G {1, 2, . . . , d + 1} and s G {+, — } in such 
a way that the number of selected boxes is minimal. See Section |4] for a more 
detailed discussion of this issue. The algorithm for calculating deg (id, [— 1, l] 2 , 0), 
I = I and s = 1 is displayed in Figure |H 

If the input of the Deg function is a sign list representing a boundary of 
an oriented cubical set, then the list of selected boxes is exactly the set Ai_ s 
from Lemma [TJ We will prove in Section 13.31 that the list faces can be split 
into J2j 9Dj+ where {Dj}j is the decomposition of Ai jS into oriented cubical 
sets and ti contains each box x the same number of times as —x, where — x 
represents the box x with opposite orientation. We will prove that Deg(faces) = 
Deg(J2j dDj) = J2j Deg(dDj). Together with equation ([2]) this implies 

Deg(L) = s(-l) l+1 ^ Deg(dDj) = deg (/, Q, 0). 

3 

One example of a possible faces construction is displayed in Figure [5l 
3.2 Pseudocode 

Function Main 
Input: 

B: oriented n-box 

/: interval-computable continuous R"-valued function s.t. ^ f{dB) 
Output: the degree deg (/, B, 0) 

boundary _info refineCov(/, B) 
return Deg(boundary_info) 

For a given sign list L and interval computable function / (see Definition [5]), 
for a box a £ L, if I(f)(a), we can infer a non-zero sign vector entry (see 
Definition [8]) for a. Moreover, due to interval-computability, I(f)(a) converges 




Fig. 3. Illustration of Lemma [T] The boundary of the oriented cubical set {_Bi,i?2} 
contains nine boxes and Ai z + = {a,b,c,d} are the boundary boxes with sign vector 
(p). This can be decomposed into oriented cubical sets D\ — {a,b} and D2 = {c, d}. 
The boundaries of D\ and D2 consist of the points marked as x and each of them is 
a sub-face of some box with sign vector different from (+) . 



New sign vector 




Fig. 4. Description of the recursive step for the identity function on the box [—1, l] 2 . 
For the choice I — 1 and s = +, we have one selected box AB. The degree functions 
returns Deg((+B , (+)), (-A, (-))) = ^"^"^ = 1. From the boxes BC and DA, 
only the sign information is used and the box CD is ignored. 



to the exact range of / on B. Hence, the following function can easily be im- 
plemented by starting with the 2n faces of B and recursively splitting the sign 
covering until interval evaluation /(/) computes the necessary sign information 
for it to be sufficient. 

Function refineCov 
Input: 

B: an n-box in M™ 

/: M™-valued interval computable function s.t. ^ f(dB) 
Output: 

Sufficient sign list wrt. /, covering the oriented boundary dB of B. 

Now, finally, we can compute the degree from a sufficient sign covering. 
Function Dcg 

Input: L: a d-dimensional sufficient sign list wrt. 

some function / covering the boundary dB of a box B 
Output: deg{f,B,0) 
if d = then 

return | ^ a sv ^ eL orientation(a) x sv 
else if L = {} then 

return 
else 

let 1 < I < d + 1 and s G {+, -} 

L sel «- {(a, (sv 1 , svd+i)) eL \ svi=s} 

L non <- L\ L sel 

faces «— {} 
for all a G L sel 

bound list of the oriented faces of a 

split the boxes in bound until for all b € bound, either 

- b is a subset of some element of L non , or 

- the intersection of b with any element of L non has dimension smaller than d — 1 
for every b 6 bound that is a subset of some box S in L non 

faces faces + (6, sv), where 

sv is the sign vector of S with omitted Z-th component 
return s(— Deg(faces) 



3.3 Proof of Correctness 

The algorithm first creates a sufficient sign list wrt. / : dB — > R" where B is 
the input box. This sign list is then an input for the recursive function Deg. We 
want to prove that if L is a sufficient sign list wrt. / covering the boundary dB 
of a box B, then Deg(L) returns the degree deg(f, B, 0). 

To prove this, we will analyze the Deg function body. When dealing with 
(f-dimensional sufficient sign lists, we always assume that some I e {1, . . . , d+ 1} 



and s £ {+, — } has been chosen. Let L be a sufficient sign list wrt. /. We denote 
L sel := {a £ L\ (SV a )i = s} and L non := L \ L sel the sub-list of selected and 
non-selected boxes. For each a £ L sel , the Deg function refines the boundary da 
until each x £ da is either a subset of some S £ L non or has at most (d — 2)- 
dimensional intersection with each S £ L non . In the first case, it assigns to x 
the sign vector SV$ with deleted Z-th coordinate. We denote the sub-list of all 
such x constructed from a by faces(a). The list faces constructed in the Deg 
function body satisfies 



and s (— Deg(faces) is returned. 

In this section, we will suppose that some implementation of the algorithm 
is given. This includes a rule for the choices of I, s, order of the lists L sel and 
L non and the choice of S. We will show that if the sign list satisfies a certain 
regularity condition defined in Definition [TUl then the Deg function output is 
invariant with respect to some changes of the input list, including any change 
of order, merging and splitting some boxes or adding and deleting a pair of 
identical boxes with opposite orientation. This is shown in Lemma [3] Further, 
we show that the list L sel can be decomposed into the sum of oriented cubical 
sets D\, . . . , D m such that ^ Uie{i m] f^l(®Dj) and such that the list faces 
constructed in the Deg function body is a merging of J2j ^Dj and a set of 
pairs {x, — x}, so that Deg(faces) = Degi^-dDj). In Theorem [3] we prove 
that Deg(J2j — J2j Deg{dDj) and combining this with equation ([2} in 

Theorem[TJ we show that if L is a sufficient sign list wrt. / covering the boundary 
dB of a box B, then Deg(L) returns the degree deg (/, B, 0). 

Definition 9. Let L and L' be two sufficient sign lists wrt. f . We say that L 
is equivalent to L' and write L ~ L' , if L' can be created from L by applying a 
finite number of the following operations: 

— Changing the order of the list, 

— Replacing some oriented box a in one list by two boxes 01,02 where 01,02 is 
the splitting of a into two oriented sub-boxes with equal sign vectors SV a = 



— Merging two oriented sub-boxes a%, 02 of a with the same sign vector to one 
list element (a,SV a ), 

— Adding or deleting a pair of oriented boxes a and —a where —a is the box a 
with opposite orientation (the sign vectors SV a and SV- a does not have to 
be necessary equal in this case), 

— Changing the sign vectors of some oriented boxes so that that sign covering 
is still sufficient and wrt. f . 

Clearly, ~ is an equivalence relation on sign lists. 




a 



SVa, = SV a 



Definition 10. Let L be a d-dimensional sufficient sign list wrt. f. We say that 
L is balanced, if each sub-face x of some a £ L such that for each b £ L, either 



iCi, or x (~1 b is at most [d — 2) -dimensional^, satisfies 

\Sx\ = \S—x\ 

where S x is a sub-list of L containing all a s.t. x is an oriented sub-face of a. 

In other words, x is a sub-face of some oriented box in L the same number 
of times as —x. 

A sign list representing the oriented boundary dB of an n-box B is clearly 
balanced, because for each (n — 2)-dimensional sub-face of some a G dB that 
is small enough to have either low-dimensional of full intersection with each 
b G dB, a is an oriented sub-face of exactly one a £ dB and —a is an oriented 
sub-face of exactly one a' G dB. The following Lemma says that the property 
of being balanced is also preserved in the faces construction procedure. This 
implies that all input lists L within the recursive Deg function are balanced. 

Lemma 2. Let L be a sufficient sign list wrt. f that is balanced. Then the list 
faces(L) created in the Deg function body is also balanced. 

The proof of this is technical and we postpone it to the appendix. 

Lemma 3. Let L be a balanced sufficient sign list wrt. f and L' be equivalent 
to L. Then Deg(L) = Deg(L'). 

Proof. We prove this by induction on the dimension of the sign list. If L is a 
O-dimensional sign list, then nontrivial merging and splitting of a box is im- 
possible. Independence of order of the list follows from the formula Deg(L) = 
\ SaeL orientation(a) x SV a and adding a pair (x, —x) to the list will add to 
the sum \ {SV X - SV. W ) = ±(sign (/ (x)) - sign (/(a;))) = 0. 

Assume that the Lemma holds up to dimension d— 1. Let L' be a permutation 
(i.e. the same multiset with different order of elements) of a d-dimensional sign 
list L and faces' be the list created for L 1 in the Deg function body. Changing the 
order of the list possibly changes the order of L sel and L non . However, a G L sel 
if and only if a G (L') sel and the same number of times. Further, faces(a) and 
faces' (a) can be constructed from each other by a finite number of splitting, 
merging and sign vector changing, because both are sufficient sign list wrt. /-,; 
representing a sign covering of the set 

U {|x| | a; is a boundary sub— face of a and x C n for some n G L non }. 

So, faces' ~ faces and Deg L = Deg faces = Deg faces' — Deg(L'). 

Further, let L' be created from L by splitting or merging some oriented box. 
If we split or merge a non-selected box, then faces'(L') will be equivalent to 
faces{L), because the equivalence class of faces(a) depends only on the union of 
all non-selected boxes. Splitting a selected box a into oi, a% will result in splitting 
some elements of faces (a), possibly changing their sign- vector (depending on the 



3 Here x and b represent just the box, without taking care of the orientation. 



choice of S in the algorithm) wrt. /-,; and generate a finite number of new pairs e 
and — e s.t. e G faces'(ai) and — e G faces'^)- So, faces(L) is again equivalent 
to faces' (L') and we can apply the induction. 

Assume that we change the sign vector of an element in L in such a way 
that we still have a sufficient sign list wrt. /. If we change the sign vector of a 
box such that we don't change a selected box to a non-selected or vice versa, 
then this change may only result in a possible change of sign vectors in faces 
wrt. f-,i (and possibly splitting and merging of the boxes in faces, if the sign 
vector change has an impact on the choice of <S* G L" on in the algorithm). So, 
in this case, faces(L) ~ faces'(L'). Assume that we change the sign vector 
SV a so that some a G L non will become selected. We claim that on the level of 
faces, this is equivalent to adding a sign list da covering a boundary of a to 
faces, all x G da endowed with the old sign vectors SV a with Z'th entry deleted. 
Choose a decomposition da — a sel U a non , where the oriented boxes in a sel have 
at most (d — 2)-dimensional intersection with each non-selected box and each 
oriented box in a non is a subset of some non-selected box (after changing the 
sign vector of a). Each x G a non was contained in some non-selected box even 
before changing a from non-selected to selected, so x G a non will be contained 
in faces'(a) list after the change. Due to the fact that L is balanced, for each 
x G a sel , there exist the same number of boxes u in L s.t. x is an oriented 
sub-face of u as boxes v s.t. — x is an oriented sub-face of v, a being among 
the u's. All the u and v are in L sel , a changed from being in L non to L sel . 
The list faces is equivalent to a list containing —x one more time then x, so 
after deleting a finite number of of pairs x, —x, faces(L) is equivalent to a list 
containing — x once. In f aces' (L'), this disappears which is equivalent to adding 
one x to the list. So, faces' {a selected) ~ faces(a non-selected) + da. We will 
show in the next paragraph that adding the full boundary da endowed with 
constant sign vectors to a sign list does not change the Deg function output. 
Summarizing this, Deg L'{a selected) = Deg faces' {a selected) = Deg faces(a 
non-selected+da) = Deg faces(a non-selected) = Deg L(a non-selected). 

Now we need to show that adding the full boundary da of a endowed with a 
constant sign vector does not change the Deg output. In the O-dimensional case, 
this says that Deg(L + da) = Deg(L) + i (s — s) where s is the sign of / on a. 
Let L' = L + da be a sign list of positive dimension such that all elements in da 
are endowed with the same sign-vector. In the consequential faces construction, 
either all boundary faces of a will be selected or they will be all non-selected. 
In the first case, faces'(L') will be a sum of faces(L) and pairs (x, — x). In 
the second case, da may be refined so that each element is either a subset of 
some other non-selected box, or has only lower-dimensional intersection with 
each non-selected box (except a) . Those a G da that are a subset of some other 
non-selected box can only possibly change the sign vector of some boxes in faces 
and boxes j3 G da that have only lower-dimensional intersection with each non- 
selected box (except a) will cause (after possible merging and splitting faces) 
the addition of a sum of pairs x and — x to faces due to the fact that faces is 
a balanced sign list. So, faces'(L + da) ~ faces(L). 



Finally, adding a pair of two selected boxes a and —a will create additional 
pairs x and —x in faces. Adding a pair of two non-selected boxes a and —a may 
enlarge the union of the non-selected boxes. Let x be any boundary sub-face of 
some selected box such that x has at most (d — 2)-dimensional intersection with 
each b G L non and x C a. Then the balancedness of the list L says that x is 
a sub-face of some selected box the same number of times as —x. So, the new 
faces are equivalent to adding each such x the same number of times as — x to 
the faces list. If we add two boxes a and —a such that a is selected and —a non- 
selected, we may change the sign vector of —a (due to the previous paragraph) 
so that both a and —a are selected and the Deg function output doesn't change. 
□ 

Theorem 2. Let B be an n-box. Then the algorithm terminates for any contin- 
uous interval- computable function f : B — > R n such that ^ f{dB) and returns 
the degree deg (/, B, 0). 

Proof. The theorem is a consequence of statement [2] of the following: 

1. Let fi\,...Qk be oriented cubical sets of dimension d + 1, let Li, . . . , Lk be 
d-dimensional sufficient sign lists wrt. to a function / : U\f2i\ — > M. d+1 s.t. 
the boxes in Li are d-boxes forming an oriented boundary dfli of fli for all 
i. Then Deg^ L % ) = ^ Deg(L z ). 

2. Let fl be a (d + l)-dimensional oriented cubical set and let L be a d- 
dimensional sufficient sign list wrt. a function / : \fl\ — > M. d+1 , such that 
the boxes in L form an oriented boundary of fl. Then Deg(L) returns the 
number deg (/, fl, 0). 

We prove both statements Q] and [2] by induction on the dimension d. If the 
sign lists are 0-dimensional, then Deg(L) = i ^2 aeL orientation(a) x SV a where 
SV a is the 1-dimensional sign vector of a 6 L and Deg(^2 i Li) = J^. Deg(Li) is 
true for any sign lists Li. For statement [21 the Deg function result is compatible 
with the one-dimensional formula 

deg (/, a%, 0) = - (sign f(b) - sign /(a)) 

for / : at ->■ M. 

Assume that the dimension is d > and both [1] and [2] hold for lower- 
dimensional sign lists. First we prove [3] Let L be a sufficient sign list such that 
its oriented boxes form the boundary dfl of a d + 1-dimensional oriented cubical 
set fl. We know that L is balanced. Let I £ {1, . . . , d + 1} and s G {+, — } be 
chosen in the Deg function body. For each box a 6 L, choose an index 1(a) s.t. 

— if (SV a )i = s, then 1(a) = I and s(a) — s 

— if (SV a )i ^ s, then choose 1(a) and s(a) so that the sign vector (S'V a )/( a ) = 
s(a) ^0 

Such index 1(a) and sign s(a) exists for each a, because the sign list is sufficient. 
For each /' G {1, . . . , d + 1} and s' G {+, — }, denote Ai>, s > a list of all boxes in 



a G L such that 1(a) = I' and s(a) = s' The list of selected boxes L sel is formed 
exactly by the boxes in Ai tS and the list of non-selected boxes is L non := L\L sel . 
For each j and each a G Dj s , let faces(a) sp iit be a splitting of faces(a) such that 
for each e 6 faces(a) sp ut and each Z> £ <9J?, either e C 6 or efl& is at most (d — 2)- 
dimensional. Further, define <9-D/ s to be the sub-list of $Z og£ >4 fo.ces(a) sp iit 
containing all x such that there exists a unique a £ Dj s s.t. a; is a sub-face of a. 
This is a sign list covering an oriented boundary of Dj (see Def. [5]). 

Define faces sp u t := J2aeL^' faces(a) spHt . Let x G faces sp i it \J2j 9Dj s an d 
assume that a; G faces(S) sp ut for S € Dj s . Because a; ^ 9Dj s , a; is a sub-box of 

exactly two boxes S and S' in Z?^ and x C b for some non-selected box b. The 
sub-list faces(S') S put contains a box y s.t. y C\ x is (d — l)-dimensional (y is a 
sub-box of some face e of <S" and e n b is (d — l)-dimensional). The orientation 
of y induced from S' is different from the orientation of x (see Def. @|. So, after 
possible further splitting of the list faces sp ut, we may assume that y = —x and 
that for each j, X) a en i f aces ( a )s P ut contains either both x and —x or none of 
them. It follows that the list faces sp nt contains x the same number of times as 
—x and the list faces is equivalent to ■ dDj s . Now we derive 



Deg(L) = s (~l) l+1 Deg(faces) = (Lemma© = s (-l) l+1 Deg(^ dDj s ) = 

3 

(induction, □) = s(-l) l+1 ^ De 9( dD i, s ) = (induction, 2.) = 

3 

= s(-l) l+1 J2deg(f^,Dj s ,0) = s(-l) l+1 ]T deg(f^Dj, s „0) = 

j j; V=l and s' — s 

(Theorem H equation ©) = dcg (/, [2,0). 



It remains to prove [TJ Assume that statement [T] holds up to dimension d — 1, 
and [2] holds up to dimension <i. Let L = L^, be a d-dimcnsional sufficient 
sign lists wrt. / such that the boxes in Lj form oriented boundaries dfli of 
oriented cubical sets Qi for i = 1, . . . , k. 

In the same way as before, we define for i = 1, . . . , k D(i)j, , to be oriented 
cubical sets such that Li is the disjoint sum Ylj i> s ' s" tne s ^S n vec tors 

have Z'th component s' on D(i)j, s , and the oriented boundaries dD(i)j s are 
sub-lists of a splitting of faces(L) such that each a; G dD(i)j g is a sub-face of 
some b 6 D(i)j, , for some Z' ^ Z. Similarly as before, faces is a equivalent to 



Deg(L) = s Deg(faces) = (Lemma© = s Deg^ dD(i)j s ) = 

(Induction, QJ) = s{-l) l+1 ^ Deg(dD(i){ s ) = (Induction, [1) 

i,3 

= s(-l) l+1 J2deg(f^D(i){ s ,0)= S (-l) l + 1 £ d eff (/^, D{t)^ s „ 0) = 

i,j i,_7; and s'=s 

(Equation ©) = deg (/, Qi, 0) = (statement 1.) = Deg(Li) 

i i 

which completes the proof. □ 



4 Experimental Results 



We tested a simple implementation of the algorithm on several linear and poly- 
nomial functions / and boxes B. All timings were measured running version 1.0 
on PC with Intel Core i3 2.13 GHz CPU and 4GB RAM. Interval computations 
were done based on the library smath [11] implementing intervals with floating 
point endpoints and conservative rounding. In theory it could happen that 64 
bit floating point representation does not suffice for computing sufficient sign 
covering of dB, but in our experiments we did not find a single example where 
this happened. 

Unfortunately, up to the best of our knowledge, all published articles on 
general degree computation algorithms, only contain examples of low dimension, 
for which our algorithm tends to terminate with a correct result in negligible 
time. Hence, in order to show the properties and limitations of our algorithm, 
we chose different examples that scale to higher dimensions. 

The first part of the algorithm where boundary boxes are subdivided and sign 
vectors are computed takes usually about 5-50 times less then the combinatorial 
part where the degree is calculated from the list of boxes and sign vectors. 
However, if there is no solution of f(x) — on B (and the degree is zero), then 
the second part terminates immediately, because — in the simplest case — there 
are no "selected boxes" at all. 

In most cases, computation of deg (/, B, 0) such that £ /(£>) \ f(dB), ter- 
minated in reasonable time if dim B < 10. If ^ f(B), then the degree is zero 
and the algorithm terminates very fast even in much higher dimensions. 

Example 1. For the identity function on [—1, 1]", the degree computation termi- 
nates even in high dimensions. The times are given in Figure [5] 




Fig. 5. In this case, B is an oriented cubical set containing two 3-dimensional boxes and 
the Deg function input dB consists of twelve 2-dimensional boxes. Two of them are se- 
lected and form an oriented cubical set D with oriented boundary 3D — {a, b, c, d, e, /}. 
The list faces contains two more boxes, identical with opposite orientation. 




Fig. 6. Time spent on calculating deg (id, [— 1, l] n , 0) = 1. 



Example 2. We considered the function 

f _ 2 _ 2 _ 2 

ji — x 1 x 2 ... x„ 

/ 2 = 2xix 2 

Jri ^X\X n . 

This function has a single root in x = of index 2 for n even and index 
for n odd. Figure [7] shows the time consumed for calculating deg (/, B,0) for 
B = [—1, 1]" and B = [—0.001, 1]™. The computation is significantly faster for 
B = [— e, 1]™ where e > is small and the root is close to the boundary. In this 
case, the subdivision of the boundary contains only two selected boxes (both 
close to 0). For B = [— e, e] n , the calculation takes about the same time as for 
B =[-1,1]". 




Fig. 7. Time needed to calculate the degree deg (/, [-1, l] n , 0) and deg (/, [— e, 1]", 0) 
for e = 0.001. 



The following table shows the number of selected and non-selected boxes in 
the subdivision of dB for B = [—1, 1]™. 



Dim B 


Selected boxes 


Non-selected boxes 


5 


32 


800 


6 


64 


2368 


7 


128 


6528 


8 


256 


17408 


9 


512 


44032 


10 


1024 


108544 



If we chose the box to be [e, 0]™ or any other such that ^ f(B), the degree 
calculation terminates almost immediately even in dimension 1000 and more. 



We also investigated the effect of the choice of I and s in the Deg function body. 
By default, they are chosen so that the number of selected boxes is minimal. 
Numerical experiments show that the computation takes more time if the num- 
ber of selected boxes is larger. The following table shows the number of selected 
boxes for various I and s in dimension 6. 



1 s 


Nr. 


of boxes 


I s 


Nr. of boxes 


1 + 


64 




1 - 


1600 


2 + 


G4 




2- 


64 


3 + 


G4 




3- 


64 


4 + 


G4 




4- 


64 


5 + 


96 




5- 


96 


6 + 


96 




6- 


96 



Choosing the bad strategy choice I — 1 and s = — would increase the computa- 
tion time significantly. The following table shows the time comparison. 



Dimension 


Optimal choice of I and s 


Worst choice of I and s 


6 


0.09 s 


0.8 s 


7 


0.65 s 


10 s 


8 


5.6 s 


175 s 



In general, we made the observation that for a fixed choice of I and s, for some 
permutations of variables the number of selected boxes, and hence run-time, is 
much higher than for others. Hence, our strategy of choosing I and s makes the 
run-time of the algorithm much more robust. 

5 Appendix 

5.1 Proof of Lemma Q] 

Let us adopt the notation a ^ B for "a is an oriented sub-face of B " (see Def. 
[3]). Let dQ be an oriented boundary of the oriented d-dimensional cubical set 
Q and let Ay lS > = {a e dQ \ (SV a )v = s'}. dQ is a disjoint union of the sets 
Ai'^i, (l 1 , s') G {1, . . .,d} x {+, — }. For each box a G Ai, s , let da be an oriented 
boundary of {a} such that for each x € da, the following conditions are satisfied: 

— for each b € dQ, either x C b or x n b is at most (d — 3)-dimensional, and 

— for each B e Q, either x C B or x n B is at most (d — 3)-dimensional. 

Such a subdivision of the boundary of a can be constructed by splitting the 
boundary faces of a as long as some boundary face x 6 da has nontrivial (d — 2)- 
dimensional intersection with some b G dQ or some B G Q. Denote by dAi :S the 
set of all boxes x G da s.t. a G Ai tS and a; is a sub-face (not necessary oriented 
sub-face) of some b G Ayy for (I', s') ^ (I, s). Finally, for any oriented box a, let 
—a be the same box with opposite orientation. 

For each B G Q, we define an oriented boundary dB of B that satisfy the 
following properties: 



— dB contains all b G dfl s.t. b <^-> B, 

— If x G da for some a G df2 and x C B, then there exist unique 61, 62 G dB 
s.t. a; *->■ 61 and —x » 62- 

Such 9£? can be constructed by completing {b G 9J? | b ^ B} to an oriented 
boundary of B in such a way that x G U a <£dnda, x C B is contained in exactly 
two sub-faces 61 and 62 of B with compatible orientations. For each c G 9B, 
either c G <9i? or c has only lower-dimensional intersection with each element in 
8S7 (if c would have a (d — l)-dimensional intersection with d G df2 and d ^= c, 
then d would be also a sub- face of B due to the second condition of Def. [3] and 
c, d G 9B would violate the first condition of Def. [5]) . 

Let Z G {1, . . . , d} and s G {+,—}. We construct the sets Dj and dDj induc- 
tively by associating the boxes in A^ s with sets Dj. Assume that D\,dD\, . . . , 
Dk-i,dD).-\ satisfy the conditions 

— Di C Ai jS is an oriented cubical set 
-An D./ = for i ^ j 

— dDj C 9/1;. s is an oriented boundary of Dj 

for all i,j < k. Let C yl i s be an oriented cubical set such that Dk flA =fl 
for i < k. Let dDk be an oriented boundary of Dk s.t. dDk C U ae _D fc 9a (such 
a boundary exists, because 9a is subdivided fine enough). If dDk Q 9Ai tS , then 
condition |U from the Lemma is satisfied for each b G dDk and the construction 
of Dk is completed. If U^ =1 Z>, ^ yl; iS , we choose some a G yl; lS \ U^ =1 Z>j and 
defining a G -Dfc+i we start the construction of a new set Dfc+i- 

Assume that 9 A ^ 9/1; iS . Then there exists some x G dDk, x £ dAi iS . 
Because x G dDk, there exists exactly one ao G Dk such that x ao (Def. [5]). 
The condition x ^ 9/1;^ implies that the intersection of x with any b G /1;'. S ' for 
(/', s') ^ (/, s) has dimension at most d — 3. We assumed that ao G dfl, so there 
exists a unique box Bi G fl such that ao B\. Let us construct a sequence 
ao, ai, . . . , a p and a sequence B^, ... , B p G J? of oriented boxes such that the 
following conditions are satisfied for u = l,...,p: 

— x <^-> a„_i B M and a„_i G 9B M , 

— a; (-a„) c — > B„ and (-<z«) G 9B„, 

— B u and B u+1 have (d — l)-dimensional intersection with compatible orien- 
tations, 

— (— a p ) G 9i?. 

The boxes B\ and ao have been defined and x ^ ao M> B%. Suppose that x <— > 
a u _i » B u . Let (— a„) be the unique oriented box in 9B M s.t. (—a;) ^-s> (— a«). If 
(— Ou) G 9i7, then u — p and we are done. Otherwise, the intersection of (— a u ) 
with each b G 9i? is at most d — 2 dimensional and it follows from Definitions 
U] and [5] that (— a u ) is a common sub-face of two boxes B u and B u+ \ in f2 with 
compatible orientations. This means that (— a u ) B u and a u B u+ i, so 
x <-} a u =-> -B u+ i. For all u, a; =-> a u , in particular — x <-» (— a p ) and it follows 
that ao and (— a p ) have compatible orientations. We add the box (—a p ) to Dk- 



We will show that this does not violate any of the above assumptions and we 
redefine dD)~ so that it is an oriented boundary of Dk and dDk C U a eD k da. 

First we show that the sequence {(a„, B u )} u terminates, i.e. it is not periodic. 
Assume that it is periodic and that (— a u ) ^ dQ for all u > 0. Let p be the 
smallest integer such that (a p ,B p ) = (ak,Bk) for some k < p. There exists 
a unique a p -i s.t. x %>-i B p and exactly two boxes -B p _i and B p in 
Q containing sub-face, so (a p _i,S p _i) is uniquely determined by 

(a p , Bp). If k > 1, then this implies (afc_i, Bk-i) — (a p _i, B p -i), contradicting 
the assumption that p was the smallest such integer. If fc = 1, then x ^ ao = 
a p _i ^ B\ — Bp and ao is a common sub- face of two elements B p and £> p -i G i7 
which contradicts ao G <9I?. We showed that the sequence {(a u , B u )} u terminates 
and we may add (— a p ) to Dk- 

If (—a p ) £ Ai >s then (— a p ) £ Ay s i for some (l',s') ^ (I, s) and a; would be a 
(d — 2)-dimensional sub-face of (— a p ), contradicting the assumption a: dAi iS . 
So, (—dp) £ The box — a p is not in Dk yet, because x is contained in both 
— a p and ao and we assumed that x G dDk- Also, (— a p ) is not contained in 
any Di, i < fc. If (— a p ) G Z^, then ao would be added to Di before, construct- 
ing the sequence (— a p ), (— a p _i), . . . , (— ai), (— ao) where x =-> a u _i ^ _B U and 
a; ^-s> (— a u ) i? M for all it. Finally, ao = — (— ao) would be included into Di, 
contradicting our starting assumption Di n Dk = 0. So, adding (— a p ) to Dfe is 
doesn't violate any assumption of the construction. 

Each x G dDj is a sub-box of some b G for (I', s') ^ (I, s). However, the 
case (£', s') — (/, — s) is impossible, because fi cannot have sign s on \Dj \ and — s 
on x C \Dj\. So, i' ^ /. In this way, we construct the oriented cubical sets Dj 
such that UDj = Ai tS . This can be done for each I and s and the resulting sets 
{Dj}j satisfy all the requirements. □ 



5.2 Proof of Lemma [2] 

Assume that L is a balanced d-dimensional sufficient sign list wrt. /. First we 
define some additional notation. We say that an oriented (d — l)-box e is small 
wrt. L, if for each F G L, either e C F or e n F is at most (d — 2)-dimensional, 
where e and F represent the boxes, without considering the orientation. Further, 
we fix the notation a M> B for "a is an oriented sub-face of B " (with the induced 
orientation, see Def. [3]) like in the proof of Lemma[T] and the notation a G D b for 
"a is an oriented sub-box of b" (Def. [2]). Further, let us represent the list L as 
a set of pairs L ~ {(-Ea, 1), (£72,2), . . . , |£|)}, where £j is the i'th element 

ofi. 
Let 

A = ({e | 3 (£7, i) <E L e E and e is small wrt. L}, C Q ) 

be a partially ordered set. If L 7^ then A ^ 0, because each oriented sub-face e 
oi E £ L can be refined to small oriented sub-boxes wrt. L. Let .M be the set of 
maximal elements in A. These are exactly the elements that are an intersection 
of a face d of some E G L with a finite number of boxes in L s.t. the intersection is 



still (d— l)-dimensional but a nontrivial intersection with any additional element 
from L would decrease the dimension of the intersection (or the intersection of d 
with all E E L, if it is (d— l)-dimcnsional). It follows that M is finite. Moreover, 
each e G A is an oriented sub-box of a unique element e' in M. We define the 
equivalence class [e] of some e G A to be the set {g G A \ g C c e'} for the unique 
e' G M. For e E A, let S e be the subset of L containing all (E, i) G L such 
that e ^ B. If e C„ e' £ .M, then S e — S e >, so we may define the set 5[ e ] 
for the equivalence class [e]. The balance property says that for each e £ A, 
\S[ e ] \ = l^i-ejl- For each e G M 1 define the bijection P[ e ] : S[ e ] — > S[_ e ] in such a 
way that f[_ e ] = Pr~i for all e E A. 

Let I G {1, . . . , d + 1} and s E {+, — }, L sel be the subset of L containing all 
(E,i) s.t. (SVe)i = s, let L \ L sel be the set of non-selected boxes and assume 
that faces = X)_EeL se ' fo-ces(E) was created in the Deg function body. We will 
represent faces as a set of elements (e, (E, i)) such that e G faces((E,i)) was 
created as an oriented sub-face of (E, i) E L sel in the Deg function body. Note 
that for a particular (E,i) E L, e °-> E cannot be contained more the once in 
the list faces((E, i)), so faces(E, i) contains each element exactly once for each 
(E, i) E L sel and can be identified with a set. In this set representation of faces, 
we ignore the order of the list. However, the balancedness of faces we want to 
prove does not depend on the order of faces. 

Let (e, (E, i)) E faces and i 4 e be so that e is small wrt. faces (this 
means that for each (g, (E, i)) E faces, either x C g or x fl g is at most (d — 3)- 
dimcnsional) . Let T x be the subset of faces containing all (g,(E,i)) E faces 
s.t. x ^ g. We want to show that \T X \ = |T_ X |. Let x' C D x be so small that 
for each (E, i) E L, either x C E or x fl E is at most (rf — 3)-dimcnsional (such 
a sub-box exists, because it may be constructed as an intersection of x with a 
finite number of boxes from L). T x = T y holds for any oriented sub-box y of 
x, so it is sufficient to show \T X >\ = |T_ X /|. To prove this, we will construct a 
bijection R x : T x i T_ x i. 

Let (eo, {Eq, io)) E T x . This means that eo G faces((Eo, io)) for some (Eq, io) E 
L sel . Let ei be another sub-face of Eo s.t. x' C e fl ei and ei is small wrt. i 
(such ei exists because of the condition that x' is small wrt. L). The sub-faces eo 
and ei of Eg are oriented compatibly, so {—x') =— >• ei and ei G »4. i?o has up to 
equivalence only two sub- faces e , e\ E A containing x' so [ei] is uniquely deter- 
mined by e . If there exists some (F, j) E L non s.t. e 1 C F, then ei C D E Q 
for some such that (e[, (E ,io)) G faces and (e' l7 (i?o, «o)) G T-^. In that 
case, we define R x ((eo, (E ,io))) := (e' 1; (i?o,«o))- Otherwise, ei is not a subset 
of any non-selected box, and for each (F, i) E L, e\ ^ F implies (F, i) E L sel . Let 
(Ei,h) := P [ei] ((E Q ,i Q )). We know that (^i,ii) G L set and x' ^ -e x ^ E 1 . We 
again find a box e2 in i?i such that the intersection e\ C\e2 contains x' and (— ei) 
and e2 are oriented compatibly, so — x' ^ ei <^-> E\. In this way, we construct 
a sequence of boxes ej and elements (Ej,ij) such that —x' e^+i ^ £Jj for 
J > 0, — ej £Jj for j > 1, all ej are small wrt. L, P{ e .}((Ej_\, ij-i)) = (Ej,ij) 
and ej is not a subset of any non-selected box for j = l,...,p. If e p+ i is a sub- 



set of some non-selected box, then e p C e' p+1 G faces((E p ,i p )) and we define 
Rx((eo, (Eo,i )) := (e' p+1 , (E p ,i p )). 

It remains to prove that R x is correctly defined, i.e. that for some finite peN, 
e p+ i will be a subset of some non-selected box, and that R x is a bijection. First 
we show that this procedure terminates. Assume, for contradiction, that the se- 
quence {[ej], (Ej,ij)}j is infinite. Because there exists only a finite number of 
(Ej,ij) G L and only a finite number of [ej], the sequence is periodic and let k 
be the minimal index such that ([e^], (Ek, ik)) = ([ei],(Ei,ii)) for some I < k. 
If I > 0, then (E u ei) = P [e! ]((£;_i, and = P[- e! ]((^, ij)) 

due to the assumption Pj ei ] = P[Z e ,j- It follows that £?;_i is uniquely deter- 
mined by ([ei],(Ei,ii)) and = (Ek-i, ik-i)- From the construction, 
we know that —a;' e; ^ Ei-\. However, in £7-1, there exists up to equiv- 
alence a unique — e;_i £7;_i s.t. x' (— ej_i) So, we proved that 
([e;_i], (E1-1, ii-i)) = ([ek-i],(Ek-i,ik-i)), contradicting the assumption that 
k was the minimal index with such equality. It remains the case I = 0, = eo. 
However, e is a subset of some non-selected box, contradicting the assumption 
that for each i > 0, e, is not a subset of any non-selected box. 

Finally, note that if R x (e , {E , io)) = (e p +i, (E p , i p ), then R_ x (e p+ i, (E p , i p ) = 
(eo,(Eo,io)) because each ([ej],(Ej,ij)) is uniquely determined by [e^+i] and 
(Ej + i,ij + i), so starting with (e p+ i, (E p ,i p )) will just reverse the order and we 
will eventually come to some e~o s.t. e~o is a sub-face of (i?o,io), e~o is small wrt. 
faces and eo is a subset of some non-selected box, so e~o is a the oriented sub-box 
of a unique (e , (i?o,«o)) G faces. This proves that = R~ x and that i? is a 
bijection. □ 
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